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Abstract. A Dalitz plot analysis of the OZI rule violating decay Df into ir~-ir + ir + is presented using 
different partial wave approaches. Scalar and vector waves are described by 7^-matrices; their pro- 
duction is parameterized in a P-vector approach. Alternatively, Breit-Wigner amplitudes and Flattc 
parametrization are used. Special emphasis is devoted to scalar mesons. The /o(980) resonance provides 
the most significant contribution. Adding /o(1500) to the scalar wave leads to an acceptable fit while 
introduction of /o(1370) and/or /o(I7IO) does not lead to significant improvements. A scan of the scalar 
wave optimizes for M = 1452 ± 22MeV/c 2 . When /o(1710) is added, the mass uncertainty increases, 
and the fit yields M — 1470 ± 60MeV/c 2 which is fully compatible with the nominal /o(1500) mass. 
The scalar wave seems to exhibit a phase motion of 270° units in the mass range from 1200 to 1650 MeV/c 2 . 
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1 Introduction 

Decays of charmed mesons provide an efficient tool to 
study meson-meson interactions at low energies. Outstand- 
ing examples are the analyses of the D + — > 7r~7r + 7r + 
Dalitz plot which helped to establish the er(500) PQ, and 
the two reactions D° -> A^tt+tt" [2] and D+ -> K~tt + tt + 
[2114] which revealed the existence of the k(700) meson. 
The reaction — > tt~it + tt + was proven to provide ac- 
cess to mesons in which the primarily formed ss state 
converts into a ir~Tr + pair [5l[6l[7l[8] in a OZI rule violat- 
ing transition. In D and D s decays into three pseudoscalar 
mesons, a large fraction of the cross section is assigned to a 
pseudoscalar meson recoiling against a scalar meson; this 
fact makes D and D s decays very well suited for investi- 
gations of the spectrum of scalar mesons and their flavor 
wave function. A survey of data sets and of Dalitz plot 
analyses of D and D s decays [9] can be found in the Re- 
view of Particle Properties [10] . 

In spite of the large potential for illuminating contribu- 
tions to our understanding of scalar mesons up to a mass 
of ~ 1700McV/c 2 , the situation concerning the states 
/o(1370) and /o(1500) and their contributions to the re- 
action Df — > 7r~7r + 7r + is still controversial. The evidence 
for the existence of /o(1370) has often been questioned 
[TT | IT2 l [T3" l ll4|: if it exists, the flavor decomposition of the 
scalar states remains unclear as evidenced by the numer- 
ous mixing schemes in which a scalar glucball is supposed 
to intrude into the spectrum of scalar qq meson, to mix 
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with them thus creating the observed pattern of scalar 
resonances. Instead of a 'narrow' / (1370), a wide scalar 
background has been proposed which was called /o(1000) 
by Au, Morgan and Pennington [12], and 'red dragon' by 
Minkowski and Ochs [13J . A survey of different mixing 
schemes and a critical discussion of their foundations are 
reported in a recent review [14) . 

All analyses of decays into tt^tt + tt + agree on ba- 
sic features of the data even though different partial wave 
analysis techniques were applied which led - in impor- 
tant details - to rather different conclusions [5][H][7][5] . In 
all analyses the Dalitz plot is shown to be dominated 
by /o(980). There is possibly a small contribution from 
tt + p°; significant contributions stem from 7r + /o(1450) and 
7T + /2(1270). In the scalar isoscalar partial wave, there is 
a sizable contribution which leads to a peak at a mass 
of about 1450 MeV/c 2 . The origin of this enhancement 
is however controversial. E791 fitted the Df — > ir + n + ir~ 
data with a scalar resonance for which mo =1434 ± 18 ± 
9MeV/c 2 and lb = 173±32±6MeV/c 2 were found. They 
had observed /g(1370) production in an analysis of the 
D — > 37r Dalitz plot and identified the scalar intensity with 
/o(1370). The /o(1370) is not supposed to have a large ss 
component, but its production could be assigned to the 
annihilation diagram (see[T];). The Focus collaboration re- 
ported m =1475±10MeV/c 2 and A = 112±24MeV/c 2 . 
These parameters are nearly compatible with /o(1500). 
Thus, the important issue to which extend the two states 
/o(1370) and /o(1500) contribute remained unsettled. 

In this paper we report on a further study of the Df —> 
tt~tt + it + Dalitz plot. The main point of the analysis is an 
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exploration of different methods and of the impact the 
analysis technique has on the final result. 

The paper is organized as follows. In section [2] we re- 
mind the reader of some basic features of weak decays of 
charmed strange mesons. Subsequently, a short survey is 
given of the apparatus (section [3]), of the data and of data 
selection. Section |4] contains a description of the ampli- 
tudes used to fit the data, and the fit results. The paper 
ends with a discussion of the results and a short summary. 

2 Decays of charmed mesons 

Before starting a partial wave analysis of Df — ► tt^tt + tt + 
decays (inclusion of the charge conjugate reaction is un- 
derstood), it may be useful to remind of some properties 
of weak decays of charmed mesons relevant to light-quark 
spectroscopy. 

Fig- [H depicts Feynman diagrams for the decay of char- 
med mesons. In Cabibbo favored decays of D+ mesons 
(Fig. [TJt) the primary c quark converts into an s quark 
under emission of a W + while the s quark acts as a spec- 
tator particle. The most likely decay of the W + produces 




(a) (b) (c) 

Fig. 1. Decays of charmed meson: a) Leading mechanism. The 
qq pair can be ud or us. b) Color suppressed diagram. The qq 
pair is likely ud. c) Annihilation diagram. Here, q can be u, d 
or s quarks. 

7r + . Conversion of a c quark into a d quark or W + decay 
into K + are Cabibbo suppressed. Hence most likely, a me- 
son with hidden strangeness emerges, recoiling against a 
7r + . Evidence for this diagram can be found in a compar- 
ison of Df — > tt + (/) and — > tt + lu. The former reaction 
occurs with (4.4 ±0.6) ■ 10~ 2 frequency, the latter reaction 
with (3.4 ± 1.2) • 1(T 3 : the reaction D+ — > it + lu is sup- 
pressed. The suppression is not as large as expected from 
cj>-uj mixing (for a deviation Sy ~ 3.5° from the ideal mix- 
ing angle, a suppression by arctg~ 2 5y ~ 1/250), suggest- 
ing other mechanisms to contribute to Df — > t: + lu. These 
branching fractions and a few further ones are collected in 
Table □ 

Fig. [TJa shows a Cabibbo allowed but color suppressed 
diagram. The two quarks of the ud pair need to match in 
color with the s and the s quark, leading to an expected 
1/N C reduction of the probability. However, Df mesons 
are observed to decay into K + K*° and K* + K° with a 
mean rate of (7.9 ± 1.4)%, to be compared with the (4.4 ± 
0.6) • 10~ 2 frequency with which Df — > ir + cf) is produced. 
Thus additional processes contributing to production of 
open strangeness are required. 

The annihilation diagram of Fig.[TJ; contributes to Df 
decays with an a priori unknown fraction. The intcrme- 



Table 1. Selected Dt decay modes [TO] 
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diate W + may convert into the 37r final state via pn + 
or /2(1270)7r + but also into a scalar and a pseudoscalar 
meson. 

Evidence for annihilation can be deduced from leptonic 
decays modes [10]. Df mesons have a large probability 
to decay into t + v t ; decays into or e + v e are sup- 

pressed due to helicity conservation. The u, d (constituent) 
quark masses suggest hadronic branching ratios to con- 
tribute several % when the number of colors is included. 
_D+ — > Tr + cu production is thus likely due to the annihi- 
lation diagram. The strange mass is of course larger; an- 
nihilation could be the reason for the unexpectedly high 
K + K*° and K* + K° yields. Cabibbo and color suppressed 
decays could also contribute to 7r + w production by rescat- 
tering processes in the final state; these processes are how- 
ever expected to make at most small contributions. 

Large fractional yields in purely pionic final states are 
observed via production of r\ and rj in the intermediate 
state. It is clear why these yields are comparatively large: 
due to their ss components, rj and rj 1 couple strongly to 
the initial ss system. Their nh component leads to pionic 
final states. The large yields of scalar mesons in D+ decays 
points to the presence of both, nn and ss components in 
the wave functions of scalar mesons. 



3 Detector and data analysis 

3.1 The .E791 experiment 

The data analyzed here have been collected at Fermilab 
in a 500 GeV/c -k~ beam impinging on platinum and car- 
bon targets. The pion beam was tracked in proportional 
wire chambers (PWC's) and silicon microstrip detectors 
(SMD's) in front of the targets; particles emerging from a 
hadronic reaction were detected in a spectrometer which 
consisted of further PWC's and SMD's, two magnets, 35 
drift chamber (DC) planes, two gas Cerenkov counters, an 
electromagnetic calorimeter, a hadronic calorimeter and a 
muon detector composed of an iron shield and two planes 
of scintillation counters. A full description of the detector, 
of data reconstruction and of data analysis can be found 
in [![T5]. 

In a first analysis stage, events were required to have 
a reconstructed primary production vertex whose loca- 
tion coincided with one of the target foils. Furthermore, 
events had to have a well separated secondary decay ver- 
tex (more than 4<7 ve rtex rocon error in the longitudinal sep- 
aration), and 3 reconstructed tracks with a total charge 
+1. Particle identification was not required, tracks were 
assumed to originate from charged pions. The mass dis- 
tribution of events due to three-pion systems surviving 
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Fig. 2. a: The 7r _ 7r + 7r + invariant mass spectrum. The dot- 
ted line represents the sum of the expected contributions from 
D° -> A:~7r + and Df — > r;'7r + decays; the dashed line is the 
total background. Events in the hatched areas at the D s mass 
are used in this analysis. The hatched area at the D + mass 
was used to establish the <r(500) and its properties, b: Detec- 
tor acceptance in the mass range from 1.97 to 1.99 GeV/c 2 . 



these cuts is shown in Fig.[2]for 7r~7r + 7r + invariant masses 
falling into the 1.7 to 2.1 GeV/c 2 mass bin. The data can 
be described by two Gaussians at 1.87 and 1.97GeV/c 2 
with resolution 01.87 = 11.8 and 01.97 = 13.5MeV/c 2 , re- 
spectively. For the background function, an exponential 
form was chosen. There are 1172 ± 61 events due to D — > 
7r~7r + 7r + events and 848 ±44 events due Df — > 7r~7r + 7r + . 

After a cut in the tt^tt + tt + invariant mass between 
1.95 and 1.99 GeV/c 2 as indicated in Fig. Hi, 937 events 
arc retained for further analysis. The integrated signal to 
background ratio was estimated to ~ 2. 



3.3 The background 

The fit in Fig. [2^ assigned 568 events to D s decays and 
280 events to the background. The latter was studied us- 
ing Monte Carlo simulations. The most prominent back- 
ground contributions stem from 

1. a general combinatorial background, 

2. D + — > K~7T + ir + decays where a K + is wrongly inter- 
preted as a 7r + , 

3. Df -> ?/7T+ decays followed by rf -> p°(770) 7 , p°(770) 
— > 7r + 7r _ with a missing 7, 

4. D° — > K~ir + decays with an extra track (mostly from 
the primary vertex) . 

Amount and shape of the background is determined using 
both Monte Carlo simulations and data. The contribution 
of the Df — > ?/7r + background is negligible. The sum of 
contributions 1, 2, and 3 above, Bgi, populates uniformly 
the Df —> 7r~7r + 7r + Dalitz plot region. In the observed 
Dalitz plot, its distribution is proportional to the accep- 
tance. This latter background represents 88 ±2% of the to- 
tal background contribution to the Dalitz plot. The shape 
of the background from D° —> K~ir + decays, Bg 2 rep- 
resenting 12 ± 2% of the total background, is taken from 
simulations. It can be described analytically, the function 
is reproduced in the second part of the r.h.s. of cq. {!]). 
The total background, Bg = Bg\ + Bg 2 , is a function of 
the (symmetrized) Dalitz plot variables: 

Bg(s 12 ,s 13 ) = 1.103 ■ Acc{s 12 , s ls ) + (1) 

0.0065- (ffD (si2,Sl3) +5D (si3,Si 2 ))/2 , 



3.2 The detector acceptance 



where g Do (s 12} s i3 ) = 

7 2343 e -0- 5 ((0- 25 -sig.)/0.2) 2 + ((3.15-si3)/0.11) 2 ) 



The detector covers the full solid angle and, to a good 
approximation, the acceptance is flat over the Dalitz plot. 
On the other hand, some of the cuts introduce small bi- 
ases and the exact shape of the acceptance needs to be 
known. The acceptance was determined from the recon- 
struction efficiency of D s production and decay and the 
phase space which is smeared out by the finite detector res- 
olution. The reconstruction efficiency was derived by the 
E791 collaboration using a full Monte Carlo simulation, 
from the tt — N interaction to the digitalization and se- 
lection of the events. The generated events had the nomi- 
nal D s mass, the mass distribution of reconstructed events 
was compatible with the distribution shown in Fig. [2^,. For 
the analysis presented here, this distribution was divided 
into 100 slices; then events were generated, for each slice, 
simulating Df — > tt^tt + tt + decays uniformly spread over 
the phase space. The summation of all Dalitz plot mul- 
tiplied with the reconstruction efficiency gave the accep- 
tance Acc(si2, S13) as presented in Fig.[2}3 (S12 = m 2 _ + , 

TT 7T (1) 

In this way, edge bins are properly taking 



Sl3 



T (2) 



After summation over the Dalitz plot, the background 
Bgi (Bg 2 ) represents 246 (34) events. The Dalitz plot dis- 
tribution of the total background is shown in Fig. [3^. 
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Fig. 3. a: The background Dalitz plot, b: The Df -» 7r _ 7r+7r + 
Dalitz plot. Since there are two identical particles, the plot is 
symmetrized. 
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3.4 The D+ — > tt-tv+tz+ Dalitz plot 

The symmetrized Dalitz plot for the Df candidates shown 
in Fig. [3)3 is reproduced from [6]. Since the two ir + are 
identical, the Dalitz plot is symmetrized. To avoid prob- 
lems with statistical errors, all entries receive a weight \. 

The most striking features of the Dalitz plot are the 
narrow horizontal and vertical bands just below syi ~ 
S13 ~1 GeV 2 /c 4 and a complex structure at 2GeV 2 /c 4 . 
The bands correspond to /o(980)7r^ decays, with / (980) 

— > 7ri,7r _ . The two ir + are identical; two interfering am- 
plitudes contribute to the final state. The structure at 
2GeV 2 /c 4 contains contributions from /2(1270)7r + , from 
p°(1460) 7r + , and from 7r + recoiling against scalar inten- 
sity. The clarification of this structure is the prime aim of 
this paper. 



4 Partial wave analysis 
4.1 The fit function 

The partial wave analysis describes the Dalitz plot of Fig. 
[3}j by a summation over possible reaction mechanisms for 
P/+ — > 2it + it~ decays plus background contributions. The 
different reaction mechanisms can interfere, hence they 
are represented by amplitudes. We consider the following 
reactions: 

1. A 7r + recoiling against 7r + 7r~ in 5-wave 

2. A 7r + recoiling against 7r + 7r _ in P-wave 

3. A 7r + recoiling against 7r + 7r~ in P-wave 

The P-wave is always described by a /2(1270) Breit-Wigner 
amplitude. The S- and P-waves are alternatively repre- 
sented by a sum of two-channel Breit-Wigner amplitudes 
or by PJ-matrices. Isotensor interactions are not included. 

Due to Bose symmetry, the amplitudes are symmetrized 
with respect to the exchange of the two ir + mesons. We 
use the following notations: P = k\ + k 2 + k 3 is the mo- 
mentum of the Df meson (P 2 = s); ki, k-i and k 3 are the 
pion momenta of the ~k~ and of the two 7r + . Obviously, 
k\ = k\ = k\ = to 2 holds. The amplitude depends of two 
invariant energy variables: 



S12 = [k\ + k 2 ) 2 and s i3 = (ki 



(2) 



The total amplitude includes contributions of different 
partial waves (S-, P- and D- waves); it can be written 
in the form 



A t ot(si2,Si 3 ) = 

J2 a [Z(k u k 3 )A a {s 12 ) 



(3) 



Z(ki, k 2 )A a (si 3 )} 



where Z is a function referring to the spin-orbital angular 
momentum structure of the two mesons emerging from 
the decays of an intermediate state. For scalar, vector and 
tensor waves, they have the following form (i = 2, 3): 



Zs{ki, h) = 
Z T {kx,ki) = ^(3^ 



1, 



Zli 



Z v (ki,ki) =zu, (4) 
[feTPH ■ 



The amplitude contains the summation over the dif- 
ferent isobars. Each isobar depends on free fit parameters 
like masses, widths and coupling constants. 

The intensity observed in the Dalitz plot is propor- 
tional to the squared reaction amplitude, to the phase 
space density <P for D s decays and the background, and 
to the detector acceptance Acc. 



N(si2,si 3 ) = \Atot(si2, S13) | 2 • Acc(s\2, Siz)<Pd s (5) 
+a-(B gi (l-p)+(3Bg 2 )<P Bg 

The two background functions are given in eq. (P). Note 
that the background definition includes the acceptance. 

For the amplitude Atot(si2, S13) we use different approa- 
ches. The scalar amplitude contains a choice of a series of 
poles for <7(500), / (980), /o(1370), / (1500), and/ (1710). 
The vector amplitude comprises up to three poles at 770 
MeV/c 2 , 1460MeV/c 2 and 1740MeV/c 2 . The tensor wave 
is described by a Breit-Wigner amplitude for the /2(1270) 
meson. In most cases, we use fixed masses and widths. 

As dynamical function describing meson resonances, 
we use if-matrices or relativistic Breit-Wigner amplitudes. 

1. In a first approach, we use if- matrices to describe 
scalar resonances. For vector mesons, either a X-matrix 
or Breit-Wigner functions are introduced. 

2. In the second method, the /o(980) is represented by 
the Flatte formula. The total amplitude is formed as 
a sum of Flatte and Breit-Wigner amplitudes. 



4.2 Breit-Wigner formalism. 

4.2.1 Breit-Wigner amplitudes 

The use of Breit-Wigner amplitudes offers large flexibility 
and good control of the fit ingredients. The price is the 
violation of unitarity when two Breit-Wigner amplitudes 
overlap. 

In the Breit-Wigner approach, the dynamical ampli- 
tude Atot(si2, s i3) is given by a sum of relativistic Breit- 
Wigner amplitudes which can be written in the form 



Ar = BW R (s) 



P = Pr 



Ml- 
-f (-)" 

PR J \PRj 



irM R 

B' l (p, P r) 



(6) 



where Mr is the mass and Pr the partial width of the reso- 
nance R. The complex number = A D +^g n n is given by 
the product of the production amplitude A D +^ of the res- 
onance R in D s decays, and the coupling constant q-kk 

for 

its decay into tttt. The (running) decay momentum is p = 
i yj (s — 4m 2 ) while the decay momentum at the nominal 
mass of the resonance is given by pr = \ \J (APr — 4m 2 ). 
The phase space is written as p(s) = \/(s — 4m 2 )/s and 
p R {M 2 R ) = v/(M 2 - 4m 2 )/M 2 , respectively. 
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For the three lowest orbital angular momenta, the Blatt- 
Weisskopf factors B' L have the form 



B' (p,p R ) = 1 
B'Mpr) = 



(zr-3) 2 +9zr 
(z-3) 2 +9z 



1 + ZR 
1 + Z 

(7) 



where z = (\p\d) 2 is the meson radius. In the fits, the scale 
parameter d is restricted to d < 0.8 fm. 

4.2.2 The Flatte parametrisation 

The /o (980) is taken into account using the Flatte parame- 
trisation 



4.3 K-matrix approach 

4.3.1 The formalism 

As reference fit, we use a if-matrix in the P- vector ap- 
proach of Anisovich and Sarantsev [16j . applied to the 
analysis of D and D s decays in [7] . The S- wave amplitude 
has the form: 



I - iptjKij{s\2) 



(10) 



where 



BW 



/o(980) 



[s) 



X 



(7T7r) 

/o(980) 



AKK) 
l ^ A /o(980) 



Sl2 



M 2 

JW /o(980) 



s - iglpw - ig 2 K PKK 



(8) 



where M /o(980 ) and I/ o(980) arc the / (980) mass and K t3 (s 12 ) = <^ 



width, A^Tjg s and A^\£gL are complex numbers. The first 
number (A 71 ""') is the same as the one used in the Breit- 
Wigner parametrization. The second expression (X KK ) 
refers to non-resonant two-kaon production and rescat- 
tering into 7T7t. The coupling constants for / (980) decay 
into 7T7T and KK are denoted as g^ and gx', the two- 
particle phase space is written as = y (s — Am 2 )/ s 
and p Kk = y/(s - Am 2 K )/s. 

4.2.3 Parameters of resonances 

The T-matrix poles of the most important resonances are 
listed in Tables[5]and[3] In most of the fits described below, 
masses, widths and some coupling constants are frozen. 

Table 2. The p (1450) and / 2 (1270) mass and width (M + i£) 



.(KK) 



Ml 



Sl2 



^scatt 
lij 



1 GeV 2 - 4 catt 

Sl2 — Sn 



p(1450) / 2 (1270) 
1460 + i 150 1275 + i 92.5 



(si2 - s A m 2 2){l - s A0 ) 

7 ^ ^ °ij 

(Sl2 - sao) 

The production of a two-meson intermediate state j is 
enhanced due to formation of resonances. The summation 
in pop extends over all considered scalar resonances a 
with mass M a . The enhancement is proportional to their 
production strengths A a (for which we use S a , V a , and T a 
to denote the strength to produce scalar, vector or tensor 
resonances), and their couplings to the decay channel 
j. The constants dj allow for nonresonant production of 
that channel; however, this flexibility is not exploited in 
this paper. The rescattering series is summed up in the K- 
matrix. Transitions from unobserved channels i into the 
observed channel j by resonant rescattering are included, 
nonresonant meson-meson scattering processes are taken 
into account by matrix elements c^. The pij form the 
phase space matrix elements: 

Pij = -yj{s- (m, + m-j) 2 ) (s - (uii - rrij) 2 ) (11) 



Table 3. Masses and widths (M + i-j in GeV/c ) of scalar 
resonances used in fits. 

<j(500) /o(1370) /o(1500) /o(1710) 

0.55 + i 200 1.32 + i 155 1.490 + i 58.0 1.769 + i 170 



The /o(980) is described by the Flatte formula; the 
following parameters are used [16j : 



M 



/o(980) 



1.023 GeV; 



(9) 



9l 



(0.15- 0.20) GcV^; ^=2.4 
9i 



In all cases, we keep the names of Particle Data Group 
even though we assign different masses to the particles. 



The rescattering process KK — > tttt is allowed below the 
KK threshold; in this kinematical region, p^ becomes 
imaginary. The parameters /fj att and SQ Catt describe a 
slowly varying part of the X-matrix elements. The masses 
M a and the couplings g™, /j S J att are process- independent 
properties; they were determined in |16| from a large num- 
ber of data sets and enter as fixed parameters into the 
analysis presented here. The constants A a and dj deter- 
mine the dynamics of a production process; they are free 
parameters to be determined from the fits. 

In a scattering situation, Chiral Perturbation Theory 
forces the 7T7t amplitude to vanish at small energies. This 
is a kinematical effect which is not enforced in production. 
The suppression of the 7T7t scattering amplitude at small 
energies is taken into account by the Adlcr- Weinberg zero, 
(s- s A m%/2)(l- sao)/(s- sao)- We choose s A = 1, s A0 = 
-0.15. 
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For higher partial waves, only the resonant terms are 
required to describe data. The AT-matrix now contains 
Blatt-Weisskopf barrier factors Bl 



#«o»ia)=x; 



1 



1 



B' ., Ml - S12 B' ., 



which depend on the orbital angular momentum L be- 
tween the two mesons, the c.m.s. momentum of the two 
mesons in the initial and final state and one scale pa- 
rameter. The Blatt-Weisskopf factors are given in eq. j7|). 
The production vector for channel Pj is multiplied by one 
Blatt-Weisskopf factor taking into account the angular 
momentum barrier for channel j. 

In this paper, the AT-matrix from [16] is used for the S- 
wavc which included five resonances. These five resonances 
are observed when the T-matrix poles of the scattering 
amplitude are inspected. The T matrix is defined by 



T = (I - iK ■ pY x K 



(12) 



Table 4. A-matrix parameters for the scalar isoscalar wave 
using 5 A-matrix poles. Masses and coupling constants are in 
GeV. Only the i = 1 /j? terms are listed since they are the 
only values relevant to the three-pion decay. 







9kk 




9vv 




0.65100 


0.22889 


-0.55377 


0.00000 


-0.39899 


-0.34639 


1.20360 


0.94128 


0.55095 


0.00000 


0.39065 


0.31503 


1.55817 


0.36856 


0.23888 


0.55639 


0.18340 


0.18681 


1.21000 


0.33650 


0.40907 


0.85679 


0.19906 


-0.00984 


1.82206 


0.18171 


-0.17558 


-0.79658 


-0.00355 


0.22358 


scatt 


/■scatt 
111 


/■scatt 
J12 


/■scatt 
Jl3 


/■scatt 
Jli 


/■scatt 
.115 


-3.92637 0.23399 


0.15044 


-0.20545 


0.32825 


0.35412 



Table 5. A"-matrix parameters for the scalar isoscalar wave 
using 4 A-matrix poles. The pole representing /o(1370) was 
omitted resulting, in [16] . in a increased x 2 - Masses and cou- 
pling constants are in GeV. 



The T-matrix poles depend on the AT-matrix masses and 
the couplings g" and are thus independent of the pro- 
duction process. The poles were identified with /o(980), 
/o(1370), / (1500), / (1710) of the Particle Data Listings 
[TU] plus a broad resonance /o(1470) which was interpreted 
as scalar glueball. 

AT-matrix masses and coupling constants g" and /jj att 
are taken from [TB] ; the values are reproduced in Table 2] 
and [5] for two cases. In Table [4] a solution is given which 
includes a pole for /o(1370); the parameters in Table [5] 
are optimized for the case where /o(1370) is omitted. The 
AT-matrix parameters reproduce the T-matrix poles given 
in Tables [land [3] 

For the vector mesons, a A'-matrix is used which con- 
tains p(770), p(1450), and p(1770). The AT-matrix masses 
and coupling constants are reproduced in Table [6] 



4.4 Fit results 



4.4.1 Study of vector and tensor contributions 







9kk 


ff47T 


9vv 


9t}t\' 


0.69566 


0.66859 


-0.71636 


0.00000 


-0.10133 


-0.17721 


1.24458 


0.87180 


0.59038 


0.00000 


0.38115 


0.39085 


1.54280 


0.36228 


0.24167 


0.50442 


0.15035 


0.45126 


1.84274 


0.07290 


-0.10118 


-0.30778 


0.03018 


-0.17852 


^scatt 


/■scatt 
111 


/■scatt 
ill 


/■scatt 
J13 


/■scatt 
714 


/■scatt 
Jl5 


-3.84414 0.54458 


-0.10980 


-0.55817 


0.46870 


-0.02877 



Table 6. A'-matrix parameters for the vector isovector wave. 
Masses and coupling constants are in GeV, d is given in fm. 
Only the i = 1 dj terms are listed since they are the only 
values relevant to the three-pion decay. 



m a d 


<?7T7r 


9kK 


g™ 


0.77873 0.34469 


0.69461 


-0.50000 


0.00000 


1.50000 0.80000 


-0.20000 


0.54000 


-0.65181 


1.90000 0.80000 


0.18042 


-0.71729 


-0.30664 



Our fit strategy is as follows: we first describe the S- wave 
and P-waves using the ^-matrices with fixed pole struc- 
ture as given in Tables [4] and [6] The /2(1270) is included 
as Brcit-Wigncr resonance. This fit has 13 parameters to 
describe the 77 independent cells. Hence there are 64 de- 
grees freedom. The x 2 of the fit is 63.49; the fit quality 
is very satisfactory. Figs. [4] and [5] show mass square dis- 
tributions and, respectively, the x 2 per Dalitz plot cell. A 
negative sign is attached to the x 2 whenever the fit result 
exceeds data. 

Results of this fit and the x 2 achieved is given in Table 
[7] as solution 1 (Sol. 1). The (complex) production ampli- 
tudes are given without errors; the statistical errors are 
small in comparison to the spread of results when the fit 
hypothesis is varied. 



In a next step wc explore the role of vector and tensor 
mesons maintaining the full flexibility of the S'-wave de- 
scription. First, we replace the AT- matrix by three Breit- 
Wigncr amplitudes. A slightly improved x 2 is obtained, 
5x 2 = —0.85 (Sol. 2). The results do not change sig- 
nificantly when p(770) is removed from the fit. Without 
p(770), x 2 increases by 0.21 compared to solution 2, the 
number of parameters by 2 (Sol. 3). Removing p(1770) 
changes x 2 by 1.75 (Sol. 4), and removing both, p(770) and 
p(1770), leads to x 2 = 64.41 (Sol. 5), an increase in x 2 by 
about 2 units while four parameters are spared. We con- 
clude that p(1460) is sufficient to describe the contribution 
of vector mesons to the data. However, if the latter reso- 



E. Klempt, M. Matveev, and A.V. Sarantsev: Further study of decays into tv tt + -k^ 



Table 7. Fit parameters for five different solutions exploring 
the role of vector and tensor resonances. The full mr S-wave 
is always included. Si and S2 interfere to make the /o(980), 
S3 = /o(1370), S 4 = /o(1500), S 5 = /o(1710). There are up 
to three vector resonances, Vi = p(770), V2 = p(1450), and 
Vz = p(1770), and the tensor Ti = / 2 (1270). Given are the 
complex production strengths A a of eq. (|10|l and the x 2 of the 
fit. 



Sol.l 



Sol. 2 



Sol.3 



Sol. 4 



Sol. 5 



Si 

S3 

Si 

s 5 

Vi 

v 2 
v 3 



8.10-i5.03 
-7.83+i2.25 
-1.81+i5.46 
-2.22+i5.95 

4.94-i6.13 
0.10+i0.14 
3.00-U5.03 

7.98-i7.76 
-5.43+i0.00 



8.61-i3.80 
7.71+il.76 
2.34-H6.85 
1.78-H7.53 
6.82-i7.85 
-0.06-i0.12 
2.86-i0.23 
3.08+i0.61 
5.30+i0.00 



8.70-i3.68 

7.73+H.64 

2.54+i6.49 

2.55+i6.94 

6.91-i7.10 

2.84-i0.29 
2.90+i0.69 
5.31+i0.00 - 



6.68-i6.37 
-7.54+i3.36 
-2.01+i9.66 
-l.lO+ill.69 
12.70-il0.78 
0.04-i0.04 
2.22 +i0.75 

5.623 +i0.000 



6.78-i6.31 
-7.55 +i3.29 
-1.82+i9.28 
-0.75+ill.21 
11.63-il0.70 

2.15+i0.85 

-5.58+i0.00 



63.49 



02.54 



62.75 



64.29 



64.41 



nance or the /2(1270) is taken out of the fit, \ 2 increases 
by 20 or more. Both these resonances are required to get 
a good fit. A further small improvement is achieved when 
mass and width of the p and fi resonances are allowed to 
vary freely The resulting values remain compatible with 
PDG values. Both are included with fixed values for mass 
and width in all subsequent fits (see Table [2]). We note 
that some fits prefer to split p(1450) into a p(1250) and a 
p(1500) even though the statistical evidence for the split 
p(1450) is weak (8\ 2 = 4 for two more parameters). 



4.4.2 Study of the scalar partial wave 

Table [H] presents a series of fits in which amplitudes con- 
tributing to the 7T7r S-wave are explored. For convenience 
of the reader we reproduce solution 5 of Table[7|as solution 
6 in Table [8j We then add a constant N\ = describ- 
ing direct 7T7r production. This gives an improvement of 
2.67 in x 2 (Sol. 7). Direct production of a KK pair (N 2 ) 
and subsequent rescattering into 7T7r is unimportant: \ 2 
changes by 0.32 units only (Sol. 8). 

We now take out /o(1370). This is an alternative so- 
lution presented in |16j . In a fit to a large number of re- 
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Fig. 5. Reference fit (solution 1 of Table [TJ to the data 
Dt —* 7r + 7r + 7r~. The \ 2 contributions of the individual Dalitz 
plot cells are given in a grey scale and in numbers. For data 
exceeding the fit, the cells are dark and the x 2 i s plotted with 
a positive sign. When the fit exceeds data, the cells are white 
and x 2 i s plotted with a negative sign. 



Table 8. Study of the scalar wave in the A"-matrix approach. 
See caption of Table [7] 



Sol.6 



Sol. 7 



Sol.8 



S0I.9 



Sol. 10 



51 6.78-i6.31 6.50-i4.39 

5 2 -7.55+i3.23 -7.78+i2.31 

5 3 -1.82+i9.28 0.88+i4.43 

5 4 -0.75+ill.21 2.68+i4.32 

5 5 ll.63-il0.70 2.68-U1.03 
JVi -1.71+i4.12 
N 2 

V 2 2.15+i0.85 2.86-H0.29 

Ti -5.58+i0.00 -5.07+i0.00 



6.81-i4.22 9.05-il.07 8.42-i2.19 
-5.39-i0.13 -8.73+i0.53 -9.71+i0.36 
0.45-il.418 -1.79+il.49 -2.51-i0.62 
-1.18-i2.64 

-0.02 -i3.52 10.97-i3.34 4.73-i7.80 
-2.41+i6.10 0.05+i4.69 
1.30+il2.76 

2.80+i0.37 1.70+i0.20 2.58-i0.25 
-4.92+i0.00 -5.39+i0.00 -5.05+i0.00 
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The experimental mass distributions are compared 
Reference fit (solution 1 of TableO to — » ir + -K + ir~ 



actions, the omission of /o(1370) resulted in a worse \ 2 ■ 
The if-matrix parameters obtained in that fit are listed in 
TablcEJ In the fit to the D+ — > 7r + 7r + 7r~ data, x 2 changes 
by 6.44 units when solutions 6 (with parameter N\ = 0) is 
compared to solution 9, and by 4.38 when solution 7 with 
N\ ^ is compared to solution 10, see Table |8j Includ- 
ing N2 in the fit with no /o(1370) gives again a marginal 
improvement. 

In a next set of fits, we used the Flatte parametrization 
to describe the /o(980), and Breit-Wigner amplitudes for 
all other resonances. The three states Jo (980), p(1450), 
and /2(1270) are always included. For the S-wave, the 
maximal set of resonances includes, apart from / (980), 
the cr(500), / (1370), / (1500), and / (1710). All masses 
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Table 9. Fits using Breit-Wigner amplitudes and Flatte 
parametrization for /o(980). Si = u(500); the two lines for 
/o(980), S 2 and S' 2 , give the couplings A"""" and \ Kk . S 3 = 
/o(1370), S 4 = /o(1500), S s = / (1710), V 2 = p(1450), 
7i = / 3 (1270). 





Sol. 11 


Sol. 12 


Sol. 13 


Sol. 14 


Sol. 15 


Si 


0.73-10.14 










s 2 


3.57+i5.18 


0.09+i5.35 


1.99+i3.03 


1.66+i4.48 


1.45+i5.79 


S 2 


2.53+i4.94 


-4.62+i3.27 


-3.59+i0.45 


-2.78+i3.10 


-1.73+i5.93 


s 3 


0.96-i0.58 




-2.13-il.57 






s 4 


0.71-10.71 






0.38-i0.98 




s 5 


4.31-U.36 








3.76+i2.09 




2.79+11.89 


0.11+13.02 


1.28+U.96 


1.07+i2.46 


1.38+i2.94 


Ti 


-5.39+i0.00 


-5.77+10.00 


-5.89+i0.00 


-5.45+i0.00 


-5.40+i0.00 


x 2 


62.47 


110.18 


78.48 


70.56 


81.40 



Table 10. Fits using Breit-Wigner amplitudes and Flatte 
parametrization for /o(980). See caption of Table [9] 



Sol. 16 



Sol. 17 



Sol. 18 



Sol. 19 



Sol. 20 



Si 
S 2 

s 2 
s 3 
s 4 
s 5 
v 2 

T, 



0.73-i0.14 
3.57+i5.18 
2.53+14.94 
0.96-i0.58 
0.71-i0.71 
4.31-H.36 
2.79+H.89 
-5.39+i0.00 



1.39+i4.12 
1.33+12.50 
0.15+i0.67 
0.66-10.61 
2.72-i0.92 
1.82+U.83 
5.55+i0.00 



1.76+i4.10 
1.17+11.87 

0.56-i0.70 
2.15-U.35 
1.85+11.82 
5.55+i0.00 



1.77+i3.27 
2.60+i2.26 
-2.03-i0.38 

2.34+U.82 
1.97+U.89 
5.80+i0.00 



1.78+i3.76 
-2.97+i2.00 
-0.81-i0.40 
0.29-i0.75 

1.27+12.15 
-5.61+i0.00 



62.47 



65.53 



65.97 



73.20 



69.30 



and widths are fixed to the T-matrix pole positions ob- 
tained in [IE]. Fit results are collected in Table 151 and ITD1 

A fit which includes all resonances (Sol. 11) yields 
the best x 2 , 62.47 for 16 parameters. The / (980) is not 
sufficient to describe the scalar intensity: when c(500), 
/o(1370), /o(1500), and /o(1710) are removed, \ 2 increases 
dramatically to 110.18 (Sol. 12). The fit is unacceptable. 
We now test the x 2 changes when only one of the three 
high-mass resonances is included in the fit. Inclusion of 
/o(1500) gives the smallest yield and the largest x 2 gain 
(Sol. 14). The other two resonances (Sol. 13 and 15) need 
a large yield (to describe intensity at the wrong mass) 
and bring a smaller gain in \ 2 - We conclude that /o(1500) 
helps best to describe the data efficiently. 

Fit 16 reproduces the results of solution 11. Removing 
Si = cr(500) changes x 2 by 3.06 units, and saves two pa- 
rameters (Sol. 17). The x 2 change is certainly not sufficient 
to claim that ct(500) is present but some contribution can 
also not be excluded. Its integrated fractional contribution 
to the Dalitz plot is estimated to ~ 3%. In solution 18, 
the S3 = /o(1370) is removed additionally, with leads to a 
X 2 increase of 0.44 units while two parameters are spared; 
the /o(1370) does not help to improve the fit. If /o(1500) 
is removed (Sol. 19), x 2 increases by 7.67. In solution 20, 
/o(1710) is removed, and x 2 changes by 3.77 units. Again, 
/o(1500) has the largest impact on the fit quality. The best 
X 2 per degree of freedom is achieved when the scalar wave 
is described by / (980), /o(1500), and / (1710). 

Instead of using masses and widths from [T5], we may 
chose values given in the PDG listings. The overall x 2 de- 



teriorates by 2 units; the x 2 changes remain close to the 
ones shown in Tables l9l and fTOl When the /o(1370) pa- 
rameters are changed, mass and width can adopt nearly 
arbitrary values. Its parameters cannot be deduced from 
this data. A free fit to /o(1710) properties puts its mass 
to the limit of the phase space, with a marginal x 2 im- 
provement. 

The fractional contributions of the various isobars to 
the Dalitz plot are not well defined quantities. First, in- 
terference effects make it impossible to assign a fraction 
of a data to an amplitude. The second argument is more 
technical: the fit fraction depends on the model used. Ta- 
bles 17181 and [9j [10] show the spread of results when the fit 
hypothesis is varied. From solution 14 we estimate that 
nearly 70% of the intensity comes from /o(980), 8% from 
/o(1500), 8% from p(1450), and 15% from / 2 (1270). 



4.5 Mass scans 

In mass scans, the mass of one resonance is changed in 
steps, and the change of x 2 as a function of the imposed 
mass is inspected. In a well behaved situation, the x 2 dis- 
tribution exhibits a local minimum at the optimal mass. 

The resonance is produced with an arbitrary phase rel- 
ative to the production of another meson. When the mass 
of the physical meson under study and the Breit-Wigner 
test mass coincide (and their respective widths) , the Breit- 
Wigner amplitude and phase, and amplitude and phase of 
the physical resonance agree over the full mass range of the 
resonance. The overall phase difference of the resonances 
under study relative to the other mesons is determined in 
the fit; this phase is denoted here as ipo . When the test 
mass is changed, the Breit-Wigner amplitude and phase 
do no longer match the complex amplitude in the data. 
In this situation, the fitted amplitude adjusts its phase 
in a way that the true phase and the fitted phase match 
in the mass range in which the true resonance and the 
test Breit-Wigner amplitude have a large overlap. In a 
mass scan, the phase motion is thus approximately repro- 
duced by the fitted phase. This method was developed in 
|17j : a first application showed that the t?(1440) region 
hosts one un-split resonance [18]. The analytical form of 
a Breit-Wigner amplitude leads to an overall sign change 
of the amplitude when going far below the nominal mass 
to far above the resonance, independent of its coupling to 
other channels, independent of inelasticities. The method 
explores directly the phase of a Breit-Wigner amplitude. 

We test this idea by monitoring the phase of the / 2 (1270) . 
Fig. [6] shows a scan of the / 2 mass region. The minimum 
in x 2 reached at 1280 MeV (see Fig. Eb). The 'observed' 
phase resulting from the fit shows a strong variation as a 
function of the imposed / 2 mass, stronger than expected 
for a single Breit-Wigner. A phase variation is observed as 
expected but quantitatively, the expected phase variation 
is not reproduced. This may serve as a warning that the 
method is not without risk. 

Fig. [7^, shows the scan of a scalar amplitude. We start 
from solution 13 of Table [9j fix the mass of the scalar 
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a) M(GeV) b) M (GeV) 

Fig. 6. a: \ 2 °f a fit i n which the mass of the /2(1270) is 
scanned, b: The associated phase as a function of the mass. 

meson at pre-set values, and determine the % 2 as a func- 
tion of the pre-set mass. The scalar mass is varied from 
1.2 to 1.7GeV/c 2 . A deep minimum in \ 2 is observed at 
M = 1452 ± 22MeV/c 2 reproducing findings of previous 
analyses. At a mass of 1680 MeV, a weaker but still pro- 
nounced local minimum is seen which we tentatively iden- 
tify with /o(1710). The mass shift with respect to the ex- 
pected value is possibly due to the phase space limitation. 

The phase motion observed in the scan and shown 
in Fig. [TJd exhibits a surprise: in the range from 1.2 to 
1.7 GeV/c 2 it covers more than 2tt. The region must house 
more than a single resonance. This result encouraged us 
to introduce a /o(1710) with fixed parameters. The phase 
motion expected for two resonances (plus a small con- 
tribution from /o(980)) agrees reasonably well with the 
'measured' phase, see Fig. [TJd. 

As a next step, we include /o(1710) in the fit with fixed 
parameters and study the remaining scalar wave in a scan 
with a second scalar resonance. The result is shown in Fig. 
[5^i, the associated phase motion in Fig. [H]b. There is now 
one x 2 minimum which is wider than the central minimum 
seen in Fig.[7K. When / (1710) is added, the mass uncer- 
tainty increases, and the fit yields M = 1470±60MeV/c 2 , 
a value which is not incompatible with 1500 MeV. If in- 
stead of / (1710) the /o(1370) resonance is included, the 
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Fig. 7. a: \ 2 of a fit which includes /o(980), p(1450), and 
/2(1270) with fixed masses and widths and one scalar resonance 
the mass of which is scanned. The width of this resonance 
is fixed to 110 MeV. b: The scalar phase motion as derived 
from the fit (solid line). The fitted scalar phase motion (dashed 
line) assumes two additional scalar resonances, /o(1500) and 
/o(1710). 
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Fig. 8. a: X 2 of a fit which includes /o(980), p(1450), / 2 (1270), 
and /o(1710) with fixed masses and widths and one scalar res- 
onance the mass of which is scanned. The width of this res- 
onance is fixed to 110 MeV. b: The scalar phase motion as 
derived from the fit (solid line). The fitted scalar phase mo- 
tion (dashed line) assumes one additional scalar resonance, 
/o(1500). The 'observed' phase motion covers a wider range 
than a single resonance does. 



production strength of the second scalar resonance gets 
small, and the x 2 minimum is found at 1400 MeV. Ob- 
viously, the interference of two close-by resonances can 
mimic the observed structure but we discard this solution. 

The phase motions shown in Figs.O) and[8jD are sug- 
gestive but should be discussed with some reservations. 
The observed phase is certainly influenced by the data 
but the fit does not determine the local phase. The phase 
is derived from a comparison of data and a Brcit-Wigncr 
amplitude of finite (110 MeV) width; thus, the phase mo- 
tion is smeared out. Second, the phase is determined from 
fits in which the mass of the resonance under studied is 
intentionally detuned. It is hence not guaranteed that the 
fit does not explore different local minima. In order to 
test the method, we have made the scans with different 
widths of the test Breit-Wigner amplitude (with r = 60 
or 160 MeV/ c 2 ). No significant changes were found sup- 
porting the hypothesis that the observed phase motions 
reflect some physics content. 

Thus we go one step beyond, and try to determine the 
number of resonances by comparing the observed phase 
motions with Brcit-Wigncr phase motions of resonances 
with fixed masses and widths. In Fig. [9j the 'observed' 
phase motion of Fig. [7] is shown together with two fits. 
Both fits use / (1500) and /o(1710), fit (a) adds / (1370), 
fit (b) adds the wide background amplitude /o(1000) sug- 
gested by Au, Morgan, and Pennington [12], and by Min- 
kowski and Ochs [13J . Both possibilities reproduce the 'ob- 
served' phase shift reasonably well even though /o(1370) 
leads to an expected phase motion exceeding slightly the 
observed one while introducing /o(1000) leads to a perfect 
match. 

We have tried to perform a mass-independent analysis 
by choosing complex amplitudes for the scalar partial wave 
in each mass bin, and offering these to the minimization 
process. This mass independent partial wave amplitude 
was partly restricted to 2, 3, or 4 mass bins. The results 
were unstable; we failed to derive a mass independent par- 
tial wave amplitude from these studies. 
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Fig. 9. The scalar phase motion as derived from the fit 
(solid line). The fitted scalar phase motion (dashed line) as- 
sumes three additional scalar resonances. In (a), these are 
/o(1370), /o(1500), and /o(1710); in (b) /o(1000), /o(1500), 
and /o(1710), where /o(1000) represents the wide scalar back- 
ground discussed in the text. 



5 Discussion and conclusions 



One could dream of future high-statistics high-quality 
data. There is however the possibility to combine exist- 
ing data from different experiments. The reaction Dj — » 
7T~7r + 7r + was studied at Fermilab by the E687 (434 events) 
[5], E791 (848 events) [5] and the FOCUS (1475 events) 
[8] collaborations. The BaBar collaboration has extracted 
2900 events; preliminary results were reported in a PhD 
thesis at Bochum [15]. Including Belle results, a factor 10 
can be reached in statistics. The different background con- 
tributions would help to understand the systematic errors. 
We urge that enterprize should be undertaken. 
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We have studied the Dalitz plot for Df decays into ir~Tr + ir + . 
In agreement with earlier findings, we found that the largest 
fraction of the data stems from scalar isoscalar mesons de- 
caying into 7r~7r + , in particular the /o(980) meson. The 
,o(1460) and ^2(1270) contribute significantly as well. We 
have made the attempt to understand the conflicting re- 
sults on scalar mesons which were obtained from fits to the 
£)+ — > 7r~7r + 7T + Dalitz plot. Previous fits agree that the 
data require scalar intensity in the mass region above the 
/o(980). When fitted with a Breit-Wigner amplitude, an 
optimal mass between 1440 and 1475McV/c 2 was found. 
The masses quoted were not compatible with the mass of 
neither the /o(1370) - which we rather assume to have 
1320MeV/c 2 - nor with the / (1500) resonance. Thus it 
is often argued that both, / (1370) and /o(1500), might 
be required to get a good fit when nominal (PDG) masses 
and widths are imposed. 

We reproduce these results. However, when a high mass 
resonance, / (1710), is introduced, the / (1500) resonance 
alone is sufficient to yield an acceptable fit. The x 2 as a 
function of the assumed scalar mass of the /o(1370)//o(1500) 
develops a flat floor, and the la mass interval extends from 
1.41 to 1.53MeV/c 2 . Thus the observed mass is fully com- 
patible with the hypothesis that the standard /o(1500) is 
produced in the reaction. An additional contribution from 
a /o(1370) resonance is not required. If the existence of 
/o(1370) is assumed, its parameters can be chosen to agree 
with the 'narrow' /o(1370) of the PDG or with the 'wide' 
red dragon of Minkowsky and Ochs. The statistics we used 
is not sufficient to discriminate between the three alterna- 
tives 'no / (1370)', 'wide / (1370)', or 'narrow / (1370)'. 
The <r(500) leads to a marginal improvement of the fit, 
with a statistical evidence just above one standard devi- 
ation. In the if-matrix fits, the inclusion of the /o(1370) 
resonance was necessary in fits to a large body of differ- 
ent reactions. In the fits to the D s — > 37r data, /o(1370) 
does provide a notable improvement even though this data 
alone is not sufficient to claim its existence. 
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